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Abstract 



Detecting binary black holes in interferometer data requires an accurate knowledge 
of the orbital phase evolution of the system. From the point of view of data analysis one 
also needs fast algorithms to compute the templates that will employed in searching for 
black hole binaries. Recently, there has been progress on both these fronts: On the one 
hand, re-summation techniques have made it possible to accelerate the convergence of 
poorly convergent asymptotic post-Newtonian series and derive waveforms beyond the 
conventional adiabatic approximation. We now have a waveform model that extends 
beyond the inspiral regime into the plunge phase followed by the quasi-normal mode 
ringing. On the other hand, explicit Fourier domain waveforms have been derived that 
make the generation of waveforms fast enough so as not to be a burden on the com- 
putational resources required in filtering the detector data. These new developments 
should make it possible to efficiently and reliably search for black hole binaries in data 
from first interferometers. 
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1 Introduction 

One of the theoretical challenges faced in the analysis of data from gravitational wave (GW) 
detectors is the construction of an efficient and reliable set of templates to search for binary 
black holes (BH) and neutron stars (NS). Firstly, a fast algorithm is needed to compute 
the waveforms as the parameter space of compact binaries is quite large requiring several 
100,000 templates to cover an interesting range of systems, making it prohibitively expensive 
to store the templates digitally Secondly, the templates we employ in our search must be 
sufficiently accurate representations of the true GW signal from a BH binary so that we 
miss out, as a result of the mis-match, only a small fraction of all possible events. In 
this talk we will report and discuss some recent progress made on both these fronts: We 
now have algorithms to compute templates, both in the time- and frequency-domains, that 
take only one-half of the computational time to filter out the detector data through these 
templates. In the time-domain this is achieved by solving a pair of ordinary differential 
equations (ODEs) |J, rather than using parametric representations of the phasing formulas 
that require solving integrals of rational polynomials ||. In the frequency-domain we have 
derived explicit analytic expressions of the Fourier components || that constitute an accurate 
representation of the time-domain signal truncated at, or slightly prior to, the last stable 
orbit (LSO). In addition, to an explicit Fourier domain phasing formula we also have a pair 
of ODEs in the frequency-domain giving the Fourier phase as a function of frequency 

In the test mass approximation, the LSO occurs at a frequency /lso = 22O(2OM /M) Hz, 
while for comparable masses it is most likely to occur at higher frequencies (see, for example 
Ref. 0). Until recently, time-domain template waveforms were truncated at the LSO since 
we were totally unaware of how to continue the waveform beyond the last stable orbit. This 
time truncation, harmless as it may sound, does mean a specific modelling of the frequency 
content of the merger signal which may or may not be detrimental to extracting the true 
GW signal, depending on how the true signal behaves. Our experiments showed that if the 
merger signal terminates quickly, say over less than half an orbital time-scale, then there is 
no harm in using the time-truncated waveforms, or their Fourier-domain counterparts [|J, 
as search templates. However, recent theoretical progress in the Hamiltonian description 
of ||, and the effective one-body approach to solving the binary black hole dynamics has 
shed some light on how a binary may end its life after crossing the LSO ||. More precisely, 
one has now a reasonably good picture of the last three milli-seconds in the life of a stellar 
mass static black hole binary when two Schwarzschild black holes merge leading to a single 
spinning Kerr hole. It has been suggested that a good strategy to detect compact binaries is 
to use as search templates the waveforms based on the effective one-body approach. These 
waveforms comprise the inspiral and merger phases followed by the quasi-normal mode ring- 
ing of the newly formed spinning black hole. Such a template bank should be augmented 
by an extended set of search templates, that result from a slight variation of the effective 
one-body waveforms, so as to take into account the uncertainty that lies in the modelling of 
the late stages in the evolution of these systems. 

A combination of the sensitivity of initial interferometers and the expected compact 
binary merger rates in the Universe has meant that the first likely binary source, an inter- 
ferometer network is likely to observe, is the merger of stellar mass black holes || [|, The 
LIGO-VIRGO-GEO interferometer network has its peak sensitivity in the frequency interval 
150-400 Hz, which corresponds to a range 11-30 M & of the total mass of a system, whose 
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LSO, and hence the most dominant part of the inspiral signal, occurs in the least noisy 
band-width of the network. Thus, these are our best candidate sources. Fortunately, the 
number of templates needed to search for binaries in this range of masses is less than 10. 
Hence we will not increase the computational cost, nor will we worsen the statistics (false 
alarm/dismissal rates) by enhancing the number of templates in this range of masses. 



2 Post-Newtonian waveform 

The post-Newtonian (PN) equations of motion and wave generation formalisms, when ap- 
plied to a compact binary evolution, facilitate the computation of the two polarizations h + 
and h x of the gravitational wave emitted by the system. For a system in circular orbit, 
located at a distance r from the Earth, the two polarisations are given as post-Newtonian 
expansions in the invariant velocity v by the following expressions [0, || 



2Mr]v< 

Ai-4- v = 



where 

#j 0) = -(1 + cos 2 i)cos2<^, H { x } = -2cosisin2<^, ... (2) 

M = mi + rri2 and rj = mim2/M 2 are the total mass and the symmetric mass ratio of the 
system, respectively, 5M = mi — rri2, i is the inclination of the plane of the binary with 
respect to the line-of-sight and <p is the orbital phase whose post-Newtonian expansion is 
also known presently to order v 5 . For equal mass binaries, that is, 5M = and rj = 1/4, 
the PN amplitude corrections ff| n ' x , n > 0, are less important than they are for asymmetric 
binaries with SM ~ M and rj <C 1. Although these amplitude corrections can be significant 
in certain stellar mass binaries, which the initial ground-based interferometers are likely to 
observe, such as a NS-BH binary, it has been customary to employ only the dominant term 
ff|°x in search templates - the so-called restricted PN approximation |J. We shall continue 
to adopt this approximation in this paper, but let us note here that a careful search for 
asymmetric binaries must include the full PN signal including all the amplitude corrections 
fT0|| . In the restricted PN approximation the signal recorded by the detector is the linear 
combination 

h = F + h + + F x h x , (3) 



where F + and F x denote the antenna patterns of the detector [pTTH . 

The evolution of the orbital phase can be worked out by using the PN expansions of the 
binding energy per unit mass E(y) of the system and the gravitational wave flux F(v) emit- 
ted, in an energy balance equation, namely, F = —M(dE/dt). For two bodies of comparable 
masses in circular orbit one has 

2 

E{v) = -TL (i + E 2 v 2 + E 4 v 4 + 0(v 6 )) , (4) 
F(v) = —t- (l + F 2 v 2 + F 3 v 3 + iV + F 5 v 5 + 0(v 6 )) , (5) 
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where 



* = -£±2, E4 _8i-5W 

12 24 v ; 



1247 35ry „ 44711 9271r? Q5ri 2 „ /8191 535r7\ 

F 2 = -, F 3 = 4vr, F 4 = 1 H -, F 5 = - [ 1 7i. 

336 12 9072 504 18 V 672 24 J 

(7) 

The phasing of the orbit is given parametrically by the following pair of ordinary differential 
equations (ODEs): 

dip v 3 dv dv dE F{v) 



dt M dt dE dt ME'(v) 

where E'(v) = dE/dv. There are several different, but conceptually equivalent, possibilities 
to proceed from the above equations in arriving at an explicit phasing formula. A straight- 
forward approach is to simply substitute for E and F their PN expansions, Eqs. (§) and 
(H), re-expand the rational polynomial F/E' to the correct PN order and solve the above 
differential equations to yield an explicit time-domain phasing formula. This is the usual PN 
approximation that has also been called the T-approximant [0. Alternatively, given the PN 
expansions, or other representations, of the flux and energy, one can numerically solve the 
above ODEs for equal time steps and then use Eq. (||). It turns out that this latter method 
is accurate and computationally quite fast. Other equivalent, but numerically distinct, ap- 
proaches are discussed in Ref. [|J, to which we refer the interested reader for an exhaustive 
account and comparison. Among the different approaches introduced therein, we would like 
to discuss here one in some detail, viz the P-approximants @, [J. P-approximants make 
the best use of our current theoretical knowledge and combine that with a fast converging 
re-summation technique. We shall compare its performance against a recently proposed new 
class of waveforms, based on a different re-summation technique [|5] - the effective one-body 
waveforms - which is likely to serve as the most accurate set of search templates. 



3 Effective One-Body Approach 

As discussed above, in the standard "adiabatic approximation" the evolution of the orbital 
phase is constructed by combining the energy-balance equation M(dE/dt) = —F with either 
the PN approximation, or re-summed estimates, for the energy and flux as functions of the 
instantaneous circular orbital frequency. Recently, Buonnano and Damour || introduced a 
new approach which is no longer limited to the adiabatic approximation and expected to 
describe rather accurately the transition between the inspiral and the plunge, and to give 
also an estimate of the ensuing plunge signal. The approach of || is essentially, like |fj, ^J, 
a re-summation technique which consists of two main ingredients: (i) the "conservative", 
that is non-dissipative, part of the dynamics [effectively equivalent to the specification of the 
E(v) in the previous approaches] is re-summed by a an effective one-body dynamics that 
replaces the two-body dynamics, and (ii) the dissipative part of the dynamics [equivalent 
to the specification of the F(v)] is constructed by borrowing the re-summation technique 
introduced in 0. In practical terms, the time-domain signal at the detector is given by the 
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following expression in terms of the reduced time t = t/M: 

h(i)=Cvl(i)cos(2ip(i)), v u = 



'dip_y 



(9) 



where C is a constant for a given binary depending on the masses of the two stars, the 
polarisation of the wave, the antenna pattern and the distance to the source. The orbital 
phase (p(t) entering Eq. (0) is given by integrating a system of ODE's: 



dr OH 



Pr i Pip ) j 



dt dp 

dip 

di dp v 



w = ■^—(r,p r ,p ip ) 



dp v 
di 



F v (Q(r,p r ,p v )) 



(10) 

(11) 

(12) 
(13) 



The reduced Hamiltonian H (of the associated one-body problem) is given, at the 2PN 
approximation^, by 



H(r, p r ,p„ 



1 

V \ 



1 + 2/7 



\ 



A(r) (1 + St + V \ I - 1 



B(r) r z 



where A[r) = 1 1 — - 



The damping force F v is approximated by 



B(r) 



i / 1 _ &n 



A(r) 



(14) 
(15) 



F v = -FpnM 



(16) 



^ V 2 v ^ Fp„( v u)) is the flux function used in P-approximants to be discussed 



where F Pn {v w 
below. 

The system of equations (P~0|) — , allows one to describe the smooth transition which 
takes place between the inspiral and the plunge, while the adiabatic evolution becomes 
spuriously singular at the LSO, where E'{v^o) = 0. Ref.[§] advocated to continue using the 
above system of equations after the transition, to describe the waveform emitted during the 
plunge and to match around the "light ring" to a "merger" waveform, described, in the first 
approximation, by the ringing of the least-damped quasi-normal mode of a Kerr black hole 
0. At the moment this technique is the most complete. It includes, in the best available 
approximation and for non-spinning black holes, most of the correct physics of the problem, 
and leads to a specific prediction for the complete waveform (inspiral + plunge + merger) 



1 The 3PN version of H has been recently obtained in Ref. 
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Masses f,„ f melge f qn 

(10,10) 285 428 964 

(20,20) 143 217 482 
(30,30) 95 145 321 
(40,40) 72 109 241 




1000 



f(Hz) 



Figure 1: Power spectral density of the signal weighted down by the noise spectral density 
for effective one-body waveforms for different binary sources. 



emitted by coalescing binaries. Because of its completeness, we shall use it as our "fiducial 
exact" waveform in our comparison between different search templates. 



The initial data {r 0) p° r) p°) are found using 



'0 



1 + 2 V (Jz(r ) - 1) 
1 - 3 V /r 2 



LJ 2 



0, pi 



r 2 — 3r] 



1/2 



F v (u) 



C(r ,pl)(dp°Jdr ) 



where z(r) and C(r,p v ) are given by 

r 3 A 2 (r) 



z(r) 



3r 2 + 577 ' 



C(r,p v ) 



A 2 (r) 



7iH(r,Q,p v ] 



z\r\ 



Qrj/r 2 



(17) 



;is) 



The plunge waveform is terminated when the radial coordinate attains the value at the light 
ring ri r given by the solution to the equation, 



rt — 3r lr + 5?7 = 0. 



(19) 



The subsequent "merger" waveform is constructed as in Ref.||. 

In Fig. m we show the power spectrum of the effective one-body waveform weighted by 
the noise power spectral density of LIGO for four different systems. For binaries of mass less 
than about 30 M the signal exhibits mostly the inspiral part, for 3OM < M < 5OM we 
see the inspiral and plunge while for M > 5OM one sees the plunge and the quasi-normal 
mode ringing, the inspiral part being insignificant. 



4 P-approximants 

The PN series for the flux and energy functions, Eqs. (fjj) and ([D, are asymptotic series 
that converge only in the limit v — > 0. In our application, however, we will be using these 
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Figure 2: Orbital evolution of a binary consisting of two stars of equal masses is shown 
plotted at different post-Newtonian orders. We show the usual Taylor approximant-based 
phase evolution on the left and the re-summed P-approximant-based phase evolution on the 
right. 



series in a region where v ~ 1. It is no surprise, therefore, that one finds that the orbital 
evolution inferred by solving the ODEs in Eq. (§) differ greatly at different PN orders. More 
precisely, the orbital evolution inferred by using the post-Newtonian expansions of E and 
F in Eq. ([8]) are significantly different even at high PN orders when v is close to 1 or when 
the binary system is close to coalescence. This is shown in Fig. |2] where we have plotted 
<p(v) as a function of v at different PN orders. Also plotted in thick solid line is the orbital 
evolution predicted by the effective one-body approach discussed in the previous Section, 
which is, presently, our best guess on how the phase may evolve. In the limiting case of a 
test mass orbiting a Schwarzschild black hole the post-Newtonian expansion of the phase 
can be worked out up to order 0(v u ). Even at such high PN orders the expansion is poorly 



convergent |12|, This poor convergence means that the PN expansions of the phase are 



not very reliable and their use in constructing templates will result in missing up to one in 
every three events that a detector can observe with the aid of a more accurate template [0. 
More importantly, the measurement of physical parameters will be severely biased, in some 
cases by as much as 100 %, thereby making gravitational astronomy imprecise. 

The situation can be dramatically improved by using re-summation techniques, which 
work with quantities that are conceptually equivalent to the PN expansions, but numerically 
different from them. These re-summed quantities are often rapidly convergent because they 
can capture a certain physical effect which a PN expansion is never able to. For instance, 
on very general grounds it is expected that the PN expansion of the flux of GW emitted by 
a test particle inspiralling in Schwarzschild geometry must exhibit a simple-pole singularity 
at the location of the light-right @, namely v = l/y/3. A straightforward PN expansion of 
the flux will never be able to predict this pole while a rational-polynomial approximation to 
the post-Newtonian series, will by construction, yield a pole. 

It is often the case in physical problems, that such rational-polynomials, called Pade 
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approximants, capture the location of the true pole singularity. However, to be able to 
employ this technique one must work with the correct physical quantities so that, to the 
best of available knowledge, only a pole singularity occurs in this quantity, and not some 
other form of singularity, such as a branch point. To apply the Pade techniques successfully, 
therefore, we have had to begin from more 'basic' quantities. For instance, instead of working 
with the usually defined energy function E(v), which is, incidentally, asymmetric in the two 
masses, we defined the energy function e(v) 




which we know in the test mass limit possesses only a simple pole on the real line and is 
also symmetric in the two masses. Indeed, as shown in Ref. Q, the rational polynomial 
constructed from just the first two terms in the PN expansion^] ex 2 {v) reproduces the exact 
energy e(v). Hence, by using, in the case of binaries with two comparable masses, this very 
form energy function as the quantity on which to apply the Pade technique, we are more likely 
to produce the correct pole singularity. As in the case of energy, we use a new flux function 
as a starting point and use the Pade approximant of the new energy and flux functions in 
constructing the phase evolution. We call this approach of introducing new energy and flux 
functions and then applying re-summation techniques as the P-approximant approach. 

The P-approximant phase evolution is plotted on the right hand panel of Fig. |2] and 
again compared to the effective one-body approach. Clearly, the curves on the right panel 
are more rapidly convergent and more closer to the effective one-body phase evolution, than 
the Taylor approximant phase evolution on the left panel. 



5 Overlaps 

The visual comparison of the phase evolution discussed in the previous Section is a qualitative 
test of the accuracy of different approximation schemes. Eventually, our interest lies in the 
construction of template waveforms that are effectual in capturing the true signal. Therefore, 
of primary concern is the fraction of the optimal signal-to-noise ratio (SNR)f] captured by 
an approximate waveform. In maximising this fraction one varies both the extrinsic, that is 
the initial phase and the time-of-arrival, as well as the intrinsic, that is the masses of the 
component stars, template parameters. The SNR depends not only on the properties of the 
signal but also on the frequency response of the detector. More precisely, the fraction of the 
optimal SNR achievable by a given approximant O is given by 

0= max , ( A,X ^ =, (21) 

parameters ^ A) (X,X) 



where, given two functions A(t) and B(t) their scalar product (A,B) is defined as 

(A, B) = 2 T J£-A(f)B*(f) + C.C. (22) 
Jo S h (f) 

2 We follow the convention wherein the Taylor expansion of a quantity f(v) to order v n is denoted as 
/ T ». Thus, f(v) = / T » + 0(v n+1 ). 

3 Optimal SNR is that SNR which would have been achieved had we employed the functional form of the 
true signal as a template in matched filtering the detector output. 
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Figure 3: The overlap of T- and P-approximant waveforms with effective one-body wave- 
forms. 

Here A(f) denotes the Fourier transform of A(t) and A*(f) is the complex conjugate of A(f) 
and Sh(f) denotes the one-sided noise power spectral density of the detector. 

In Fig. |3] we have plotted the overlap of T- and P-approximants with the effective one- 
body waveforms as a function of the total mass, for equal mass binaries, and for a system in 
which one of the bodies is a neutron star. We have maximised the inner product only over 
extrinsic parameters in order to exhibit how faithful are the two approximants in reproducing 
the effective one-body waveform. The superiority of P-approximants over T-approximants is 
quite clearly brought about in this plot. For binaries whose merger takes place at frequencies 
beyond the sensitive bandwidth of a detector the disagreement between P-approximants and 
the effective one-body approach is not too great. However, for sources that merge in the 
frequency band of the detector, it is best to employ effective one-body waveforms as search 
templates. 

Having demonstrated that the two re-summed techniques give very similar waveforms we 
suggest that that the effective one-body waveforms be used as search templates. Indeed, a 
good strategy would be to employ a set of search templates based on, but encompassing a 
larger space than, the waveforms derived from the effective one-body approach. 
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